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In this paper we present results of the lowest eigenvalues of random Hamiltonians for both fermion 
^ , and boson systems. We show that an empirical formula of evaluating the lowest eigenvalues of 

random Hamiltonians in terms of energy centroids and widths of eigenvalues are applicable to many 
different systems (except for d boson systems). We improve the accuracy of the formula by adding 
moments higher than two. We suggest another new formula to evaluate the lowest eigenvalues 
for random matrices with large dimensions (20-5000). These empirical formulas are shown to be 
applicable not only to the evaluation of the lowest energy but also to the evaluation of excited 
energies of systems under random two-body interactions. 
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I. INTRODUCTION 



The 7 = 0+ ground state dominance for random Hamiltonians was discovered by Johnson et al. in 1998 [jy]. Many 
efforts have been devoted to understand this problem since then. See Refs. 0, [3] and references therein for details. 

Recently, Papenbrock and Weidenmueller considered fluctuations of and correlations between the J-dependent 
spectral widths 0, Q. By using such correlations they were able to approximately evaluate spin / ground state 
probabilities for six nucleons in a single-j {j = 19/2) shell under a two-body random ensemble (TBRE). Along this 
line, a very simple formula of the lowest energy of spin / states based on energy cenroid, spectral width and dimension 
of spin / states, was presented in Ref.^ by the present authors. Our formula (see Eq. (12) of Ref. [^) was shown 
to be applicable to evaluate statistically the lowest energy of spin / states, and proved to be very good in predicting 
spin I ground state probabilities calculated by using random two-body interactions. We should note that the idea of 
evaluating the lowest eigenvalue based on energy centroids and spectral widths was suggested by Ratcliff [3 , Vary et 
al. 0, and by Zuker and his collaborators However, predicted results by formulas in Refs. 0, SB, Hly] are less 
accurate than the formula in Ref. JjJ. 

Evaluation of the lowest eigenvalue is not only useful in studying regular structure of atomic nuclei in the presence 
of random interactions, but also a very common practice in many other fields. It is therefore the purpose of this 
article to revisit empirical formulas of evaluating the lowest eigenvalues under random Hamiltonians. In this paper we 
shall also suggest other empirical formulas of evaluating the lowest eigenvalues of random Hamiltonians and random 
matrices. 

This paper is organized as follows. In Sec. II we review results of evaluating the lowest eigenvalues by using energy 
centroids and width. It is found that the results of single-j shell and sd-boson systems follow the formula of Ref. Q 
(except d bosons). In Sec. II we also improve such evaluations by adding higher moments for various systems. In Sec. 
Ill we investigate other four lowest eigenvalues. In Sec. IV we concentrate on discussions of d bosons. The lowest 
eigenvalues of random matrices are discussed in the Appendix. 



II. EMPIRICAL FORMULA OF REF. [q 

In this paper we take the same notations as our earlier work [6] . For fermions in a single-j shell 

H= V2JTTGj A'^ X , (1) 

J— 0, even 

with A'^'^ = ■^[O'lcily and A'^ = —^[ajajY. Two-body matrix elements Gj's are assumed to follow the TBRE, i.e., 
they are a set of random numbers with a distribution function 

p{G.,) - ^exp(-G2 /2), J - 0, 2, 2j - 1. (2) 

V ZTT 

Matrix elements of H for spin-/ states can be expressed in terms of coefficients of fractional parentage (cfp's): 

Hip^ = {rimrH = E ^^ipi^j- (3) 

One can obtain eigenenergies of spin / states by diagonalizing Hjp^. Such definition can be easily generalized to 
many- J shells without confusion. 

In Ref. Q, a very simple formula of evaluating the lowest energy of Hj was suggested as follows, 

=^,-<i>(d,)a,, (4) 



where $(c?/) = ^/alndi -\- b, ui is the square root of the second order moment of eigenvalues of all spin / states, di is 
dimension of spin I states, and a and h were determined empirically to be 0.99 and 0.36, respectively. This formula 
was found to hold statistically. 

Let us first investigate fermion systems. We take four fermions in a. j — 15/2 shell, four fermions in a j = 21/2 
shell, four fermions in a j = 31/2 shell, five fermions in a j = 19/2 shell, six fermions in a j = 17/2 shell, and four 
fermions in a two-j {j — 7/2,5/2) shell. The procedure is the same as in Ref. d; first, $(fi/) = [Ei - E'f™^]/ai 
are calculated for 1000 runs of the TBRE, and second, we calculate the average of ^(d/) for these 1000 runs. We 
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plot (^((i/))^ as a function of Indj in Fig. 1. One can see that (<I>((i/))^ is close to linear correlation with ln((i/) for 
all these examples. Parameters a and b obtained by this procedure in each example are shown in Table I. They are 
close to those obtained in Ref. (i = 0.99, b — 0.36). The results of sd bosons are similar to those of fermions in a 
single-j shell or two-j shells, as shown in Fig. 1(b) and Table II. We also investigate results of a more complicated 
system: three valence protons and three valence neutrons in the (2si/2, 1^3/210*11/2) shell, as shown in Fig. 2. 




FIG. 1. (Color online) Factor ($((i/))^ versus In (d/) determined numerically for fermion and sd-boson systems, 
(a) fermion systems; (b) sd-boson systems. The dimension in these cases is 2-40. 



TABLE I. a and b values for fermions in a single-j shell and two-j shells corresponding to (a)-(b) of Fig. 1(a). 



(2j,n) 


(15,4) 


(21,4) 


(31,4) 


(19,5) 


(17,6) (2ji 


= 5, 2j2 = 7, n 


= 4) Rcf. f6] 


a 


1.04 ±0.03 


1.02 ±0.01 


1.00 ± 0.02 


1.04 ± 0.03 


1.07 ±0.02 


1.06 ±0.02 


0.99 


b 


0.28 ±0.04 


0.32 ±0.03 


0.34 ± 0.04 


0.32 ± 0.08 


0.20 ± 0.06 


0.30 ± 0.04 


0.36 



TABLE II. a and b values of sd-boson systems with n — 5, 6, 7, 8, 9, and 12, respectively. See Fig. 1(b). 



n 


5 


6 


7 


8 


9 


12 


a 


1.10 ±0.02 


1.09 ± 0.03 


1.13 ±0.05 


1.08 ±0.04 


1.10 ±0.05 


1.05 ±0.02 


b 


0.24 ±0.03 


0.25 ±0.05 


0.17 ±0.08 


0.26 ± 0.08 


0.22 ±0.11 


0.30 ±0.07 
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FIG. 2. (Color online) Factor ($(c?/))^ versus \n{dj) for three valence protons {Np ~ 3) and three valence neutrons 
{Nn — 3) in the shell (2si/2, 1^3/2: 0*11/2). The line is plotted by using $((i/)^ = alndj + b. 

Eq. ^ was shown to be well applicable for evaluating spin / g.s. probability (see Ref. for details). However, 
it is not good enough for evaluation of the ground state energy in a reasonable precision. It is therefore desirable to 
improve the formula of Ref. 6] . 

Towards this goal, we consider higher orders of moment of the eigenenergies to compensate the deviation from 
Gaussian distribution for eigenvalues. In this paper we consider the cubic root of the third order of central moment, 
denoted by (T3: 

(-. dj di di di dj dj di 

i=l k=l j = l I i=l j=l i=l I i=l 




Let us assume 



Ef'"" =Ei-C2<72 + C3a3, (6) 



where (72 = c/ in Eq. According to our numerical experiments (see Table III), the disagreement between predicted 
result of by using Eq. ([6]) and that by diagonalizing Hi can be reduced to about 1/2 on average, in comparison 

with our earlier formula Eq. ([¥]). Unfortunately, we are not able to obtain a simple formula for C3 . The values of C2 
such obtained are very close to the value of ^(rf/). 

TABLE III. Relative deviation of Eqs. ([4]) and ^ for four fermions in a single-j shell with j = 31/2 and n — A. 

^prod__g,oxact 



Error (e) is the TBRE average calculated by the value of 
Error B (e^) is obtained by Eq. 



Error A (e^) is obtained by Eq. ©, while 



d 


3 


4 5 


8 


10 


14 


20 


25 


29 


I 


50 


3 


39 


2 


4 


13 


12 


20 


eA 


0.40 


0.46 0.18 


0.28 


0.33 


0.30 


0.13 


0.15 


0.17 


es 


0.31 


0.17 0.098 


0.18 


0.15 


0.12 


0.091 


0.089 


0.10 
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III. OTHER LOWEST ENERGIES 



The success of evaluating the lowest energies encourages us to go further. Here we study eigenvalues of the first 
to the fourth excited states. Typical results are shown in Fig. 3, where we take the same form of Eq. fl]). One can 
confirm here that Eq. (4) is also applicable to evaluate the excited energies for both fermions and sd boson systems 
while the values of a and b are different from the ground state. The results are summarized in Fig. 4 and Table IV. 




FIG. 3. (Color online) [<I>(d/)]^ versus \n{d). (a) the first excited states for fermions. (a') the first excited states 
for sd bosons, (b) the second excited states for fermions. (b') the second excited states for sd bosons, (c) the third 
excited states for fermions. (c') the third excited states for sd bosons, (d) the fourth excited states for fermions. 
(d') the fourth excited states for sd bosons. 




FIG. 4. Comparison of average $(c?/)^ versus \n{d). In the figure, "g.s." means the ground state, "1st" means the 
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first excited state, "2nd" means tlie second excited state, "3rd" means the third excited state, and "4th" means the 
fourth excited state. 



It is worthy to note that our empirical formulas are based on calculations with finite dimensions (dimension is 
less than 10*). According to the linear correlation between (<I>((i/))^ and ln{di) shown in Fig. 4, the straight line 
corresponding to ground states and that corresponding to excited states (e.g., the 1st excited states) seem to cut 
across each other. Such an intersection does not occur, because the g.s. energy always corresponds a larger <I> value. 
In other words, the linear correlation between {^{dj))^ and ln{dj) is valid when di is smaller than 10''. For larger dj, 
the results should be further investigated. 



TABLE IV. a and b values for fermion systems and sd-boson systems with eigenvalues of the ground state and the 
first to the fourth excited states corresponding to Fig. 4, respectively. 



State 


g.s. 


1st 


2nd 


3rd 


4th 


a 


1.04 ±0.01 


1.18 ±0.01 


1.21 ±0.01 


1.23 ±0.01 


1.25 ±0.01 


b 


0.30 ±0.03 


-1.42 ±0.02 


-2.14 ±0.03 


-2.64 ±0.02 


-3.03 ±0.02 



IV. D-BOSON SYSTEMS 

In this Section we study very simple systems, d bosons, for which eigen-energies are linear combinations of two-body 
matrix elements. The two-body Hamiltonian of a d-boson system is given by 

Hd^Ho + J2 \V2rTTci[id^ X d^y X (d X d)T . (7) 

From Eqs.(2.79) and (2.82) of Rcf. fTo|, we have 

E= Eo + ^(4c2 ± 3c4)nd{nd - 1) ± ^(7co - lOca ± 804) [nd{nd ± 3) - v{v + 3)] 

+ ^(-C2±C4)[/(/±l)-6nd] . (8) 

For states with given /, the state with v being maximum value or minimum (depending on the sign of p' = 
Yq{7co — 10c2 ± 3C4)) is the lowest among all eigenvalues. For given and /, 

Ej^ = E{nd, I) - ^(7co - 10c2 + 3c4)v{v + 3) , (9) 

where E{nd,I) is the same for all spin I states. 

Our results of ^{dj) versus ln((i/) are presented in Fig. 5. There seems no simple relationship between $ and Ind/. 
Apparently, there is systematic deviation from linear correlation between [^(d/)]^ and ln{d). The value of ^{dj) 
seems to saturate when the dimension goes to infinity. 
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Fig. 5. Phenomenological factors [^(d/)]^ versus In (d) for d bosons with different boson numbers. The hne is 
plotted by assuming [^(d/)]^ = 0.99 In (d) + 0.36. (a) na = 9, (b) na = 18, (c) Ud = 36, (d) = 72, (e) = 108, (f) 
= 180. 



The deviation from hnear correlation between $(d/) and Ind/ in Fig. 5 originates from the distribution of eigenen- 
ergies of d bosons. The eigenvalues of other more complicated systems in this paper exhibit a Gaussian distribution 
[ist . while that of d bosons is close to a triangular (/ > n) or trapezoidal (/ < n) distribution, as shown in Fig. 6. 
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Fig. 6. Relative distribution of eigenvalues of d bosons, (a) n = 180,/ — 240; (b) n — 180,/ — 100; (c) n ~ 
360, / = 460; (d) n = 360,/ = 260. 

V. SUMMARY AND DISCUSSION 

To summarize, in this paper we have studied lowest eigenvalues of random Hamiltonians. First, we demonstrate 
that our semi-empirical formula suggested in Ref. = Ej — ^{dj)ai with = alnc? + b works well for 

various examples (so? bosons, fermions in a single-j or many-j shells) with dimension ranging from 2 to 5000 (except 
for d boson systems for which there are systematic deviations) . We also improve our formula by adding the third-order 
central moment. 

Second, we investigate eigen-energies of the excited states including the first, second, third and fourth lowest 
energies. We find that the same formula with different parameters describes statistically very well to eigen energies. 

Third, we study d-boson systems and discuss why there are systematic deviation from our statistical formula of the 
lowest eigenvalues. 

In this paper we also study the lowest energies of random matrices in Appendix. We see that there is systematic 
deviation from Wigner's semi-circle prediction [16\ when the dimension is not very large: when the dimension is less 
than 100, the lowest eigenvalue can be evaluated statistically by the same formula of Ref. @. We suggest another 
statistical formula of lowest eigenvalues of pure random matrices. 
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Appendix A The lowest eigenvalues of random matrices 

All matrix elements in this Appendix are given by Gaussian distributed random numbers. We change the dimension 
of matrices with dimension ranging from 2 to 3000. In Fig. 8(a), the dimension changes from 2 to 15. In this case, we 
see a similar result as in main text. There exists linear correlation between $^ and hid (see Fig. 8(a) ), with a and b 
are very close to those in Sec. II. This means that Eq. Q applies to general cases (statistically) when dimension of 
matrix is not very large. 

When the dimension changes from 15 to 3000, we find that there is systematic deviation from linear correlation 
between $^ and \nd. We empirically obtain 

-$(d) =a\nd/d + b: a = 1.59 ± 0.03, & = -2.00 ± 0.003. (10) 

as shown in Fig. 8(b). We see that Eq. (jTU]) is more applicable than [$((i)]^ — a\n (d) + b when dimension of matrix 
becomes large. We note that the parameter b in Eq. pO)) saturates at 2 in the large d limit. This saturation value is 
the predicted value of Wigner's semi-circle theorem. 
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Fig. 7. (Color online) Phenomenological factors ^{d) versus In (d)) for random matrices, (a) matrices with smaller 
dimension. Suppose = aln (d) + b. We obtain a = 0.97758 ± 0.02426, b = 0.37201 ± 0.04754 by fitting. These 
values are very close to those in Ref.[6l]; (b) matrices with large dimension (d = 15 — 3000). One can see a nice linear 
correlation between -$(d) and Ind/d. We obtain that a = 1.589 ± 0.028, b = -2.001 ± 0.003 in Eq. 



